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Abstract. We calculate target fragmentation in pp — > nX and 7p — > nX reactions in the meson cloud 
picture of the nucleon. The pp — + nX reaction is used to fix the pmr + form factor for three different 
models. We take into account the possible destruction of the residual neutron by the projectile. Using 
the form factor from the hadronic reaction we calculate photoproduction and small XBj electroproduction 
of forward neutrons at HERA. Here the qq dipoles in the photon can rescatter on the residual neutron. 
In photoproduction we observe slightly less absorption than in the hadronic reaction. For deep inelastic 
events (Q 2 > 10 GeV 2 ) screening is weaker but still present at large Q 2 . The signature for this absorptive 
rescattering is a shift of the da/dE„ distribution to higher neutron energies for photofragmentation. 

PACS. 13.60.Hb Total and inclusive cross sections - 11. 80. La Multiple scattering 



1 Introduction 

In the last two decades the study of total inclusive deep 
inelastic scattering (DIS) processes has allowed to extract 
important information on the structure of hadrons. Parton 
distributions have been determined and scaling violations 
have been tested to a high level of accuracy. The QCD 
improved parton model has been shown to be very reliable 
in the presence of a hard scale. 

On the other hand semi-inclusive reactions with elec- 
tromagnetic probes are still less explored. They can im- 
prove our knowledge on the inner hadronic structure. 
Through the study of new observables characterizing these 
processes we may ask more detailed questions about the 
proton. In this case the application of perturbative QCD 
(pQCD) has been restricted to high p t events, and experi- 
ments have been analyzed in terms of parton distributions 
and fragmentation functions. 

With the HERA collider target fragmentation can be 
studied in a much cleaner way than with fixed target ex- 
periments. New interest has been triggered by measure- 
ments on leading neutron production performed by the 
ZEUS and HI collaborations at the electron-proton col- 
lider These data are currently analyzed in terms of 
hadronic degrees of freedom, i.e. studying the virtual pion 
flux in the nucleon Thus one hopes to extract infor- 
mation about the pion structure function at very small x, 
not reachable in Drcll-Yan experiments. 
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In the late fifties Chew and Low already suggested the 
idea of using pions from the pion cloud of the proton as 
targets to get information on the interaction of different 
projectiles with pions ||. 

Novel theoretical tools to include forward leading parti- 
cle production in the framework of standard pQCD have 
been developed by Trentadue, Veneziano and Graudenz 
0| , who introduced a new set of nonperturbative distribu- 
tions (fracture functions) which allow to absorb collinear 
singularities at leading order in the QCD coupling con- 
stant. 

If we look at 7*p reactions in the cm-system, we can 
consider the incoming photon as a qq state || that inter- 
acts with a proton made up of two color neutral compo- 
nents, one of which is the final state neutron. This involves 
quite different aspects of the nucleon than those inves- 
tigated in deep inelastic inclusive scattering. Especially 
long range properties of the baryon can be analyzed in 
such a process. At which length does the string connect- 
ing a quark to the residual diquark break and produce 
two colour neutral objects, a meson and a nucleon? This 
question is especially important for nuclear physics, as one 
wants to know what amount of the nucleon-nucleon inter- 
action is describable in terms of meson exchange forces in 
the nucleus. 

From the high-energy point of view the determination 
of the pion structure function is central. The common in- 
terpretation of the 7*p — > nX experiment relies on the 
application of the meson cloud model of the nucleon P,[lC|, 
pd[ together with the factorization hypothesis that allows 
to separate the reaction into two steps: the fragmentation 
process and the interaction |l^,|l3). It assumes that the 
fragmentation process is universal, i.e. independent of the 
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projectile which initiates fragmentation. In the following 
we will examine the validity of this factorization hypoth- 
esis carefully. 

Let us consider the generic reaction with an incoming 
projectile a leading to neutron production i.e. ap — > nX 
(see fig. Q). In the one-pion-exchange model the differential 
cross section is given by the product of the pion flux factor 
times the total air cross section: 



da a P 



dzdpt 



Fnw(z,Pt)cr 



(1) 



The flux factor Fm V (z,p t ) gives the probability for the 
splitting of a proton into a pion-neutron system. It de- 
pends on the longitudinal momentum fraction z carried 
by the detected neutron and its transverse momentum^] 
p t . The total cross section ofJJ.(s') is a function of s' the 
center of mass sub-energy squared in the a 7r + interaction. 




ments. We will investigate the relevance of absorptive cor- 
rections in detail in order to understand the one-pion- 
exchange mechanism and the extraction of the pion struc- 
ture function. We apply high-energy Glauber theory to 
calculate and compare the screening corrections in lead- 
ing neutron production for pp and 7W p reactions. We 
especially search for differences between photoproduction 
and deep inelastic scattering in neutron fragmentation. 
To this end we follow the (^-evolution of the fracture 
functions from a kinematic region where the perturbative 
evolution of is not yet applicable to large Q 2 . We 
find similarly to ref. |l7j that absorptive corrections for 
highly virtual photons do not vanish, as indicated by data 
[p]]l^1. However this effect is smaller than what has been 
found for proton-proton interactions. This dependence on 
the projectile can be a source of factorization breaking. 
With respect to ref. jL7|, we extend this analysis to real 
photons (Q 2 — 0) and we find that the size of rescatter- 
ing corrections is comparable to that of hadronic reactions 
and bigger than that of virtual photons. 

The outline of the paper is as follows. In section 2 we 
describe the meson cloud model, in section 3 we calcu- 
late the absorptive corrections to pp — > nX. Section 4 is 
devoted to virtual and real photoproduction of forward 
neutrons. Section 5 closes with a summary and a discus- 



I 7T + 



2 Meson Cloud Model 



P 



n 



In this section we review briefly the main features of the 
meson cloud model (MCM). 

In the MCM a proton is viewed as a bare proton sur- 
rounded by a virtual meson cloud, 



Fig. 1. Picture of semi-inclusive neutron production accord- 
ing to the one-pion-exchange model. 



The usual procedure ||Q is to fix the parameters en- 
tering the flux factor from the data of leading neutron 
production in proton-proton (a = p) collisions. Then as- 
suming universality, one applies the same equation (Q) to 
virtual photon scattering (0 = 7*) in DIS. Measuring the 
differential cross section "f*p — > nX one extracts the pion 

structure function, since (riot is proportional to F£ + . 
Obviously the factorization hypothesis plays a crucial role 
in this method. 

Indeed early works 



|14J,[15,[T6| on absorption effects for 
pion-exchange mechanism in hadronic inclusive reactions 
indicate a suppression factor of the order of 50-70%. Re- 
cently a first and detailed study of absorptive corrections 
in the Regge formalism (TtJ has appeared. Many of our 
conclusions agree with the ones reached in ref. Wf\ . These 
absorptive effects depend on the projectile and are a source 
of factorization breaking. They can be comparable or even 
more important than other background contributions, al- 
ready estimated, reducing the accuracy of F% measure- 

1 We use pt = |p t |. 



IPt) 



s {bo T> 

+ J2J2f dzd ^ Pt)\B, M; z, Pt ) } 

W BM 

VS { |po T) + J2 I dzd2 Pt ( z ' P*) ( 2 ) 



b,7T°; Z,Pt 



n,n T ; z,p t ) 



where <pg M {z, p t ) is the probability amplitude to find, in- 
side a proton with spin up, a baryon B with longitudi- 
nal momentum fraction z, transverse momentum p t and 
hclicity A and a meson M, with longitudinal momentum 
fraction 1 — z, transverse momentum — pt and hclicity A'. 
We restrict ourselves to the first contributions of this ex- 
pansion in terms of Fock states. V^S is the renormalization 
constant, which is fixed by (p\p) — 1 and gives the ampli- 
tude for the bare proton. 

In the light-cone approach the amplitudes 4>ntt, for a 
proton with spin +1/2, read |10 



<t>N* 0,Pt) 



Vigo 



1 



m N (z - 1) 
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J N 



1/2,0/ \ v3go 



Pi e' 



, (3) 



"TV 



where Mjy-^ is the invariant mass of the pion-nucleon sys- 
tem, given by 



ML 



m 2 N +p 2 t ml+p 2 
z 1 — z 



rriN and are the nucleon and the pion masses; go is 
the bare pion-nucleon coupling constant and if is the az- 
imuthal angle in the transverse plane. Including the renor- 
malization factor y/S from eq. (|2|) we get the renormalized 
effective coupling g = V~Sgo [jllj, which can be extracted 
from low-energy data: we use g^/in = 13.75. 

Because of the extended nature of the hadrons in- 
volved, the interaction amplitudes in eq. (^) have to be 
modified by including a phenomenological nNN form fac- 
tor. It is important to stress here that while the vertex is 
derived from an effective meson-nucleon Lagrangian, the 
form factor is introduced ad hoc. In order to parametrize 
the form factor we need to introduce the momentum trans- 
fer t which can be expressed as follows 

t = (PN ~p' N ) 2 = ~\ [Pt + (1 - z) 2 m 2 N ] 

= (l-z)(m 2 N -M%„)+ml. 

Different models and parametrizations are available in 
the literature. In the following we compare the results ob- 
tained using the light-cone approach and the covariant ap- 
proach, the last one with inclusion of reggeization. Besides 
these two form factors we consider a irNN form factor ex- 
tracted from Skyrme-type models Jl^,^o). This form fac- 
tor leaves low momentum transfers essentially unaffected 
while suppressing the high momentum region strongly. 
The three models are: 



light-cone irNN form factor 
G(z,p t ) = exp[Rf c (r 



= exp[i?^_ m 2)/(i_ z) ] ( 4 ) 

— covariant ttNN form factor 

G{z iPt )=exp[I%(t-m*)] (5) 

— Skyrme model irNN form factor 

G(z,p t )=exp[R 2 s t]g(t) (6) 

where g(t) is a rational function of t given in the ap- 
pendix. 

The amplitudes (f> xx must be changed according to </> AA — * 
(j) XX G{z,pt) . The flux in eq. ([!]), for proton fragmentation 
into a neutron can then be calculated as 

F n „(z, Pt ) = lnJ2\0NUz,Pt)\ 2 \G(z,p t )\ 2 , (7) 
6 XX' 

where 2/3 is the isospin factor and the azimuthal an- 
gle in the transverse plane has been integrated out. The 



reggeization of the pion (relevant for z — > 1) is included 
in the covariant and Skyrme approaches by the further 
change 

^XX' ^0AA' ( ' 1 _ 2 )-a.(t) ) 

where a n (t) — a 7I (Q) + a'^t is the pion Regge-trajectory, 
with 0^(0) = and a' v ~ lGeV~ 2 . The light-cone form 
factor contains the decrease of the cross section for z — > 1 
already in the exponential. Its form, however, is a crude 
approximation and therefore we do not expect the light- 
cone form factor to be adequate for extremely large z. 

In the following section we consider pp collisions in 
order to fix the parameters i?; c and R c appearing in the 
previous equations^. 



3 Estimate of absorptive corrections in 

pp — > nX 

We consider target fragmentation reactions as stripping 
reactions in the cm-system where the projectile proton 
strips a 7r + from the target proton leaving behind a neu- 
tron. The projectile proton smashes the pion into pieces, 
while the neutron remains intact as a spectator. Any addi- 
tional interactions like A— production or p-exchange may 
spoil this simple picture and reduce the accuracy of the 
determination of the itpn vertex. This has been studied 
and the amount of such a background is estimated to be 
around 20% HQ. We will neglect these processes in our 
calculation, but we model them rescaling the one-pion- 
exchange cross section by 1.2. 

We are fully aware that this rescaling procedure is a 
poor simplification of the background contributions. A re- 
cent work [ pl| shows that an improved treatment of the 
background is important to get the right neutron distribu- 
tion at very high z values. Nevertheless we will concentrate 
mainly on the region of pt =0 and 0.7 < z < 0.9, where 
our simple rescaling gives reasonable results. 

The invariant differential cross section for the one- 
pion-exchange mechanism is (in light-cone approach^]) 



E,, 



d 3 a 



da 



d 3 p n 7r dzdp 2 

2g 2 



16tt 2 z(1-z) {M* v 



(8) 

^A\G(z,p t )\ 2 aZ- 



rn 



2 \2 
NJ 



This picture is reliable when the pion and the neutron in 
the target proton are well separated, i.e. at large z and 
large impact parameter. For small impact parameters and 
intermediate z values it must be extended to allow the 
scattering of the projectile on the neutron and the conse- 
quent screening effect (see fig. ||). 

To be more precise we distinguish in our notation be- 
tween the target (p(r)) and projectile {p(p)) proton: 



P{P)P(T) 



nX . 



For the Skyrme form factor we use the value fitted in JL9[ : 
Rl = 0.031/m£. 

3 Reggeized covariant expression can be obtained multiply- 
ing eq. (g) by (l — z )~ 2a ^W anc j taking G(z,p t ) from eq. (|). 
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Fig. 2. Picture of the collision for two different impact parameters. On the left we show a peripheral collision (large b) where 
the pion is stripped and the neutron acts as a spectator. On the right we show a more central collision where the projectile 
proton can destroy the neutron through rescattering. 



We employ the high-energy Glauber approximation to mul- 
tiple scattering which has been used for target fragmen- 
tation in heavy ion collisions in ref. p^ j. 

A more formal derivation of absorptive corrections can 
be found in fl7|| , where Regge calculus and generalized 
AGK cutting rules |^3| are applied. Our approach can be 
traced back to that of and viewed as a simplification 
of it in the eikonal approximation. 

We treat the target proton as a pion-neutron system 
(4>o ) undergoing a transition to an excited state (4> a )- The 
cross section for this process can be expressed as 



,) = J d 2 ba ^ a (b) 



= d 2 b \(<f> a [l-S^Sn\<t>o)f » (9) 



CT(0O 



where b is the impact parameter and S ab are the inter- 
action operators (see below). We assume that the proton 
state <f>o can be factorized into a system of a pion and a 
neutron 

I0o) = \ipo P ^ono) = \ipl P )\ir )\n Q ) , 

with -0g P the spatial component, specified to keep spatial 
and intrinsic degrees of freedom separated. Similarly the 
excited state with an undisturbated neutron has the form: 



where we have introduced an extra index (j) to take into 
account all possible spatial configurations for this excited 
state; \ir a r) is an arbitrary state with the same quantum 
numbers as the pion. To get the total cross section we 
sum now over all spatial configurations and over all \~K a ') 
states. We apply the closure relation and exclude the elas- 
tic contribution ( | 7Tq ) ) 



<Kb)=5>0-aO»)= £E| Jd 

n j 



y n d 3 y 7r 



lpj(yn,yir)lpo(yn, y 7 r)5'^o( b - S Tr) S 00 ( b ~ S ' 

E / dZ ynd s y^\ii}o{y n ,yn)\ 2 

X S% Q (b - 8 v )Sl P J(b - sM(b - S„)| 2 

d 3 y n d 3 y^\ipo(yn,y^) 1 2 
x [l-\S%(b-s v )\ 2 ]\S$(b-s n )\ 2 , (10) 

where y n = (s„,z„), y w = (s„-,^), s„. and s„ are the 
coordinates of the pion and the neutron in the impact pa- 
rameter plane; z n and z n are their longitudinal momentum 
fractions. Also 

1 - \Soo I 2 = i 1 1 r p *\ 2 ~ 2Rer p7r 

il-2Rer pn , (11) 



\s&\ 2 = \i-i*" 



where the profile functions r describe the respective two- 
body scatterings and are related to the scattering ampli- 
tudes in momentum space by Fourier transformation 



/w'( q ) = /rf 2 be iq b r«(b). 

2tt J 



(12) 



Let us now consider the density distribution \ipo(ym j/tt ) | 2 ■ 
We parametrize this density starting from the probability 
density to find a pion and a neutron at a certain trans- 
verse separation b re j = s„ — s x , imposing the center-of- 
mass constraint in b-space and the longitudinal momen- 
tum conservation: b cm = z n s n + z^s^ = 0, z n + z n = 1, 

\4>o(yn,y-K)\ 2 = p„7r(2„,b re ;)(5 2 (z„S n + Z w S Tr )5(z„ + Z 7r -l) . 

(13) 

Inserting eqs. (y) and d||) into eq. © and carrying out 
the 3-dimensional integration over y^, we get for the semi- 
inclusive cross section 



da pp - 



*nX 



dz 



= d 2 b 



da{b) 

dz 



with 



(14) 



U. D'Alesio, H.J. Pirner: Target Fragmentation in pp, ep and yp Collisions at High Energies 



5 



rfcr(b) 

dz 



d 2 s n - 



(l-z) 2 

x 2i?e.T p7r (b - s n ) [1 - 2Rer pn (b - s„)] 

/ ^reiPmzi^Z , fore/ 

)2i?er p7r (b + 
x [1 - 2i?erP"(b - (1 - z)b rei )] , (15) 



where we have restored the previous notation, z„ = z and 
replaced the pion and neutron coordinates by 

t>rei = 

Z 



l-z 



3.1 Pion-neutron density and profile functions in 

b-space 

In order to extract quantitative information from eq. ( p"5| ) 
we evaluate the probability density to find a pion and neu- 
tron at a certain distance inside the parent proton. The 
main idea is to start with the amplitudes in momentum 
space for the splitting of a proton into a pion-neutron sys- 
tem (see eqs. (0)) and then calculate the Fourier transform 
in two dimensions with respect to the transverse momen- 
tum. 

We obtain the amplitudes in b-space (we keep z fixed) 

where by (j) l nv we mean y/2jl(j)^l. From this we obtain 
the probability density to find a neutron and a pion re- 
spectively with longitudinal momentum z and 1 — z and 
relative transverse separation b re i : 



Putt \%i ^rel ) 



\^(z,b rel )\ 2 



In fig. ^ we show the behaviour of these densities at dif- 
ferent values of z for the light-cone and the covariant ap- 
proach. 

Concerning the profile functions we start from scatter- 
ing amplitudes with gaussian shapes, 



/ a6 (q) 



iPcxn n b 

^ttexp 



2A* 



and by Fourier transformation we get 
r ab (b) = ±-ag t A 2 ab exp 



(16) 



(17) 



where we have considered purely imaginary amplitudes as 
we are interested in the high-energy regime. 



3.2 Cross section for pp — » nX 

The gaussian dependence of the profile functions allows us 
to perform the b-integration in eq. ( |l4|) analytically, we 



have then 



daPP- 



*nX 



dz 



J ' d 2 b re ip n7r (z, b re i) Cr£^ + 



(18) 



{l-A 2 



eff^exp 



A 2 h 2 



}■ 



with 



A 2 



A 2 A 2 

pn pn 



A 2 

P7T 



A 2 

pn 



\A 



A 2 pn w 0.08 GeV 2 



(19) 



Before giving explicit results from eq. Jig), let us try to 
interprete the physical picture emerging from it. 
The first term alone is nothing else than the standard ex- 
pression according to the factorization hypothesis, eq. (Q), 
i.e. the stripping of the pion cloud inside the target proton. 
The second term represents the screening correction which 
is the most interesting result of our calculation. The Born 
fragmentation cross section is multiplied by the proba- 
bility that the projectile proton does not destroy the neu- 
tron component of the target proton. The screening factor 
changes the simple factorization picture. 

Assuming that the final state interaction does not mod- 
ify the transverse momentum distribution of the fragments, 
we can calculate the invariant differential cross section 
E n by multiplying the differential cross sections for 
the longitudinal distributions with the transverse prob- 
ability distribution for the fragmentation process. This 
method seems to work quite well in nuclear target frag- 
mentation p2|. We get 



d 3 crPP- 



>nX 



d 3 p n 



z 1 dN(z, Pt ) da p P- 
7T N(z) dpi 



iX 



dz 



(20) 



where da/dz is given by eq. ([18]). The normalized fraction 
of fragmentation processes in the interval p t , p t + dp t is 
obtained from the normalized pion flux factor: 



1 dN(z,p t ) F nn (z,p t ) 



N(z) 



dpi 



Jdp 2 F n7I (z,p t ) 



(21) 



Definitely the experimental cross sections for p t — 
[[m],[2^] have the most pronounced shape. Here one can 
really see a fragmentation peak which seems to be su- 
perimposed on some background. We estimate additional 
contributions^] to have a similar shape as the calculated 
cross section and lead to a 20% correction. Thus we scale 
our result for the one-pion-exchange contribution by a fac- 
tor 1.2. A model calculation of background processes can 
be found in ref. pl[ . For larger p t values the fragmen- 
tation cross sections become rapidly flatter and smaller 
and the background increases. Therefore we fit the ra- 
dius parameters to the pt — data. In fig. || we show the 
semi-inclusive cross section for pp — > nX as a function of 
the longitudinal momentum fraction z of the neutron to- 
gether with the experimental data from ref. f24j]. For the 



4 These may come from resonance excitations like p ■ 
and p — > tt + N*°, which decay into neutrons. 



7T+A° 



6 



U. D'Alesio, H.J. Pirner: Target Fragmentation in pp, ep and yp Collisions at High Energies 




total cross sections cr%* t and o\™ t we adopt the fits^ per- 
formed in ref. fl27|. We adjust the radius parameters Ri c 
in the light-cone and R c in the covariant form factor to 
the data and find a reasonable agreement with the fol- 
lowing values: Rf c = 0.2 GeV~ 2 and R\ = 0.05 GeV -2 . 
An important feature of the screening correction is that 
its inclusion reduces the radius parameters so extracted. 
The shape of the Skyrme form factor does not deviate ap- 
preciably from the covariant form factor, because even at 
smaller values of z the momentum transfers \t\ in fig. ^are 
not large. One should also keep in mind that the reggeizcd 
pion has a variable spin different from zero, when t is dif- 
ferent from the pole value. The coupling of virtual pions to 
the nucleon does not have to be identical to that imployed 
in one-pion-exchange potentials. 
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Fig. 4. Invariant differential cross sections for neutron pro- 
duction at p t — calculated according to eq. (20|) for three 
different form factors. Data points are from ref. 



5 In fig. |E| for a^ t we take directly the scattering data at 
small energies from the particle data group fed. 



Fig. 5. Invariant differential cross sections for neutron pro- 
duction at p t — 0.2 GeV calculated according to eq. ( |20| 
for two different form factors. Skyrme-type approach, not dis- 
played, gives almost the same result as the covariant one. 
Points are low-energy y/s — 6.84 GeV bubble chamber data, 
ref. |25l. 



The differences of the light-cone and covariant ap- 
proaches are more pronounced for very large z and at 
Pt 7^ . Here the light-cone form factor vanishes faster 
with z compared to the fall off in the reggeized covari- 
ant form factor. In fig. |^ we compare the light-cone and 
the covariant-Regge models with the invariant differential 
cross sections measured at pt — 0.2 GeV [^f|. The theory 
does not agree very well with these data. In fact these 
data are taken at y/s — 6.84 GeV, where Feynman scal- 
ing is not expected to hold. One sees for the low-energy 
data that near z = 0.9 the ir + p resonances affect the cross 
section. There are also data at larger energies |2£|, but 
these do not cover simultaneously the large z and small pt 
range. 

The same calculation performed for the pp —* nX re- 
action can be applied to pn — > pX, where the same prnr 
vertex factor enters. The only changes to be done are cr p7r 
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6. Inclusive cross section for the pn — > pX process at 
= 195 GeV for \t\ < 1.4 GeV 2 , see text. Data points are 

a. 



for a p7r+ and a pp for <r pn in eq. ( |l~8| ) . In fig|j| we show the p 2 - 
integrated cross section for the pn — > pX in light-cone and 
covariant-Regge pictures with the data from p9| . Again a 
reasonable agreement is found with the values of the radius 
parameters fitted on the leading neutron data at pt = 

(fig!- 

The effect of screening is shown clearly in fig. where 
we plot the if-factor, i.e. the ratio of the differential cross 
section with and without absorptive corrections for the 
fitted values of the radii. The individual models differ 
slightly in the z < 0.8 region where there is a sizeable, 
> 30%, screening. The light-cone form factor leads to a 
slightly smaller K-i actor. For very large z values large 
transverse separations between the neutron and the pion 
dominate [O (see fig. ||), thus the projectile proton misses 
the neutron for almost any finite impact parameter and 
the .K'-factor approaches unity. 
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Fig. 7. Absorptive iGfactors for neutron production in 
proton-proton collisions for three different form factors. 



A comparison of the data with a calculation neglecting 
the screening mechanism would give larger radii. Such a 
fit can represent the data down to lower z values like z = 
0.65. In general the amount of absorption can be reduced 
if one uses a larger cut-off radius in the pion form factor 
and additional exchanges fl2l[ . An understanding of the 
trade-off between these two effects necessitates that they 
are considered together. 



4 Cross sections for 7*p — > nX and "yp — > nX 

In this section we consider the interaction of an initial vir- 
tual or real photon with a proton leading to neutron tar- 
get fragments. Before entering the detailed calculation of 
the semi-inclusive reaction it is worthwhile to make some 
remarks on inclusive photon-nucleon interaction in con- 
text with our irNN form factor. Recently p(^] an esti- 
mate of the antiquark distributions in the nucleon has 
been performed with the help of the Sullivan formula, 
cf. eq. (Q). The calculation shows that for the Skyrme form 
factor used above the MCM prediction exhausts or only 
slightly exceeds the light sea quark distributions u and d 
in the proton at Q 2 — 1.2 GeV 2 . This is a large improve- 
ment with respect to harder form factors often used be- 
fore, which strongly overestimate this contribution. Thus 
the problems discussed e.g. in ref. [jllj seem to be well 
taken care of. Two comments are in place: perhaps the 
F2 calculations are not valid because of the interference 
between photon-nucleon and photon-pion inelastic inter- 
actions, where the target-like slow fragments of the struck 
sub-hadron and the residual hadron interact |l7j] . Since 
both are slow and strongly interacting there is no rea- 
son that this final state interaction is negligible. A sec- 
ond remark about the pion pole calculation is the possible 
contribution of diffractive excitation of the proton feeding 
the neutron channel. This problem has been investigated 
in pof . Preliminary experiments [fl8| rather point to the 
fact that for z < 0.8 there is a sizable contribution of 
7To-exchange to forward proton production |3l|] . 

A new approach j2^] has been recently applied to de- 
scribe in a consistent and unified framework proton frag- 
mentation and the flavour asymmetry in the proton sea. 
This work |2l) revises old form factor parametrizations 
and background contributions to one-pion-exchange in com- 
parison with previous formulations in order to reproduce 
the new E866 data on the u — d asymmetry. 
We are aware that the models for the form factors we use 
here follow strongly the old approaches and therefore can 
not reproduce the new E866 data as well. As we limit 
our calculations to z-values bigger than 0.7, the issue of a 
quantitative description of the u — d asymmetry is not so 
relevant |3l| . 

In the semi-inclusive reaction we limit ourselves mainly 
to forward neutrons in photon induced reactions at small 
XBj, as they are studied at HERA. In this kinematic region 
one cannot apply the usual factorization into a photon- 
quark cross section times a distribution function multi- 
plied by a quark fragmentation function. That is why it 
is interesting by itself to study target fragmentation as 
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a nonperturbative process combining fragmentation and 
structure information. With good reason the new concept 
of fracture functions has been invented for this process. 
For small ccej we can consider the photon as a quark- 
antiquark state which materializes long before it reaches 
the proton || and interacts with the pion and neutron 
in the proton wave function. Schematically we write the 
inelastic cross section in a form similar to the proton in- 
duced one: 



where we have introduced the notation of average to get 
simpler expressions in the following. The second term con- 
tains the correlated average of qqn and qqn cross sections. 
Equation (E5L) reads 



tf> a ) = J d 2 bdwd 2 r\V q q{w,r)\ 2 



(22) 



Here the qq pair wave function is represented by x P q q{w, r) 
with w being the momentum fraction of the quark and r 
the transverse separation of the quarks. This wave func- 
tion contains the summation over the electric charges of 
the different quarks and is electromagnetic in origin. The 
weak electromagnetic coupling determines the cross sec- 
tion of the photon with the pion. 

Following the same procedure for the sum over all spatial 
configurations and over all excited (inelastic) states we get 



da ~<*p^nX 

dz 



where we take into account the dependence on the qq-sae 
by an r-dependent cross section: 



r'''''Mb.r) = i-a?r(r)4 g - a exp 



h 2 A 2 - 



(24) 



By performing the b-integration analytically we get 



dz 



= J d 2 b re ip njI (z,b re i) J dwd 2 r\& qq (w,r)\ 2 



x a t 



I- A 



2 o ffi(r) 
2tt 



exp 



A 2 h 2 n -i 

yl cff U rel 1 

2 if 



with A e g = 



A 2 - A 2 - 

qq-K qqn 



A 2 - 



A 2 - 

qqn 



(25) 



(26) 



In this form we see the role played by the photon wave 
function: It governs the integral but does not enter in 
the magnitude of the screening correction with the same 
weight as in the direct term. Screening is a strong inter- 
action effect which is a function of the transverse size of 
the qq pair. 



In the expression (25) we have a linear and a quadratic 



term in the dipole-hadron cross section <7** t (r) averaged 
on the photon wave function squared. The first is easily 
calculated giving the total photon-pion cross section, 



■ : : ' / — 



drfr|^>,r)|V t T(r) 



7 7T 



(27) 



da 7 P" 



dz 



d 2 b r 



:lPn 



(28) 



Ccff 



1-^cff^exp 



A 2 h 2 
/l cS °re 



where we have defined 



Ccff 



(29) 



It is worthwhile to remind here that the total cross sec- 
tions appearing in eqs. ( p8| ) and ( p9| ) depend on Q 2 , the 
photon virtuality, and on the scaling variables x v and x n 
(respectively for a pion and a neutron target): 



1-2 



^Bj 

z 



with iBj 



Q 2 



2q-p 



(30) 



= J d 2 b rel p n7I (z,b re i) J d 2 bdwd 2 r\& qq (w,r)\ 2 

x 2iteT** r (b-s Tl r)[l - 2iter 9 « n (b - s„, r)] , 

(23) 



Aware of this we rewrite eq. ( p9[) in the following way 



qqn 



Ccff 



qq-K 



qqn 



(•£•71 



qqn 



(x n )) 



Ott)) WtotM) 

(31) 

dc of'™, keeping the right x- and Q 2 - 
dependences. The quantity in eq. ( |3l| ) can be then be 
parametrized following Kopeliovich and Povh PH as 



where we use of^ 



(o-toT(^)) 



N 



1 



{%7T ) \ •^■■K 



1 



^off 



1 



( 32 ) 

This parametrization is valid in the region of large Q 2 
and small x values^. We choose the values No = 2 GeV~ 2 
and A e g — 0.15 slightly different from those quoted in 
|H as we are interested in a region of smaller x values. 
The low value for A c g comes from the fact that diffraction 
and therefore also shadowing are dominated by processes 
softer than those dominating F% for which we use the 
parametrization in the double leading log approximation 
[[34j which covers a wide xbj and Q 2 range. 

In the case of real photons at very high energies (HERA 
kinematics) we neglet the small effect in the difference be- 
tween the photon-pion and photon-neutron cm-energies 
and get 



0"efif|<2 2 =O 



(KT) 2 > 



16tt d<j~>P^ pX 



jp^X 



a 



dt 



t=o 



16tt 



20 mb . 



(33) 



where we use da lv ^ vX jdt = &£>ct£> exp(bot) with the ex- 
perimental values for bo and an according to fl35|, |36| . 

6 Working at small we can exchange all neutron labels 
against proton labels. 
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Fig. 8. Comparison of the energy distributions of neutrons 
for deep inelastic events and photoproduction imposing the 
kinematical cuts of ref. Q. The distributions are normalized 
to the total number of DIS events. 



virtuality of the photon than e.g. the free dipole-neutron 
cross section. Still screening is reduced for high Q 2 . To 
show this we plot in fig. |^ the integrated neutron energy 
distributions in photoproduction with Q 2 < 0.02 GeV 2 , 
5 GeV < E e < 22 GeV, and in deep inelastic scattering 
with 10 GeV 2 < Q 2 < 100 GeV 2 , 0.04 < y < 0.95, the 
ZEUS cuts (0 sca t < 0.6 mrad, pt < 0.5 GeV and an inte- 
grated luminosity of 6.7 pb _1 ). A shift of the DIS cross 
section maximum by 50 GeV to lower energies is clearly 
visible. This is due to screening, which reduces the cross 
section mainly at smaller E n (neutron energy), such that 
the peak appears at higher energies in photoproduction. 
Our results indicate a clear signal for more transparency 
of the neutron when interacting with highly virtual pho- 
tons. Small size color dipoles rescatter less on the target 
fragment neutron, even if for very large Q 2 screening still 
persists. The semi-inclusive reaction with Regge-exchange 
also allows to study the same phenomena as in the diffrac- 
tive cr 7 p^p+ x cross section. 



Finally following the same idea of shadowing as a soft 
process which motivated the choice of A c g we employ 
slope parameters A%g n and A 2 S7r which are calculated from 
total cross sections]] a 
20 mb: 

A 2 ~ 4 " 



qqn 
tot 



*qqn 



qqn 



A 2 - 

qqn 



30 mb and 

47T 



qqn 



For the double differential j*p cross section, we obtain 
with the help of eqs. (20) and (pl|), 



d 2 al ' P ^nX 1 dN(z,p t ) da^'P-* 71 * 



dzdp 



N(z) dp 2 



dz 



(34) 



which is then rescaled to yield the differential cross sec- 
tion for neutron production in electron-proton scattering 
(besides a small longitudinal contribution) 

d i ep^e'nX ^2 7 *P-nX 
- 2E — (2 - 2y + 2y 2 ) 



dxBjdQ 2 dzdpf 2ttxb)Q 



dzdpj 



(35) 

with y = jJtj! s is the cm ep total energy squared. 

The phenomenological j*tt + cross section entering eq.(p8[) 
is expressed in terms of the pion structure function 



7 7T^ 



47r 2 a 

~~Q 2 ~ 



(36) 



In the following we use for demonstration purpose 
F£ + (x Bh Q 2 ) = ^F%(xbj,Q 2 ) valid at XBj -C 0.1. For 
real photons we refer to the fit performed in . 

Averaging the colour dipole-neutron cross section to- 
gether with the dipole-pion one, the decrease of the ef- 
fective cross section a e s is less steep with increasing the 



7 From the neutron target to the pion target we scale total 
cross sections with the ratio of pion to proton size, i.e. with a 
factor 2/3. 
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Fig. 9. Absorptive A"-factor for neutron production in real 
and virtual photon-proton collisions. 



The absorptive J-C-factors for jp — > nX and ep — > e'nX 
are shown in fig. O: it is visible that for DIS neutron pro- 
duction (high Q 7 ) screening becomes very weak , while 
for real photons we get, at large z, a somewhat reduced 
effect compared to proton-proton collisions (see fig. 0). 
For z > 0.75 and Q 2 > 10 GeV 2 the factorization of the 
cross section into a pion flux factor and a pion structure 
functions looks very acceptable. 

In fig. [To] we show the pt-integrated neutron cross sec- 
tion with the HI cuts as function of z. The screening in 
pp — > nX manifests itself as a bigger ep — > e'nX cross 
section, since a higher pion flux is needed to explain the 
measured hadronic cross section. Notwithstanding the va- 
lidity of the impuls approximation in high-Q 2 processes, 
factorization is broken in soft hadronic fragmentation re- 
actions. The comparison of our result with the preliminary 
HI data is very encouraging. 

The fracture functions ||(7| introduced to describe tar- 
get fragmentation allow a model independent evolution in 
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Fig. 10. Comparison of the differential cross section including 
screening corrections in light-cone approach with preliminary 
HI data, ref. 0]. 



the framework of perturbative QCD, albeit the starting 
functions have to be known. In this approach the non- 
perturbative distributions are related to the pt-integrated 
differential cross section for leading neutron production in 
ep scattering in the following way 



d 3 a £ P- 



*e nX 



2-na 2 



dxBjdQ 2 dz XBjQ' 



l -(2-2y + 2y 2 )M 2 (x Bh Q 2 ,z), (37) 



where at leading order in a s there is no contribution from 
the fragmentation of the struck quark in the target region. 
In the one-pion-exchange model without screening correc- 
tions one can interpreted] M 2 (xb], Q 2 , z) as the product 
of the flux of neutrons integrated over p t , times the pion 
structure function F 2 (x-^^Q 2 ). 

The fracture function M 2 is defined in analogy with the 
structure function F 2l 



M 2 (x,Q 2 ,z) = x^2e 2 M hn/p (x,Q 2 



(38) 



where M i n / p (x, Q 2 , z) represents the probability of find- 
ing a parton of flavour i with momentum fraction x and 
a neutron with momentum fraction z inside a proton. 

In fig. [ll] we give the fracture functions M 2 (x B j, Q 2 , z) 
as functions of z, at fixed XBj, for various Q 2 values. 

The (Revolution of the fracture functions we use is 
given at small virtualities by the higher twist effect of 
screening and at higher Q 2 by the evolution of the pion 
structure function. In fact for neutron fragmentation the 
Altarelli-Parisi evolution of the fracture function dom- 
inates the total Q 2 -dependence as has been argued in 
ref. §7ll. 



Here one should consider this fracture function as an input 
distribution at a certain scale Qq and then apply the pQCD- 
evolution. 



5 Discussion of fragmentation results and 
validity of the factorization hypothesis 

Finally we compare the screening corrections for the three 
different cases of proton, real and DIS photon induced 
semi- inclusive fragmentation reactions (cf. figs. 0, One 
sees that both proton induced and real photon induced 
cross sections have if-factors differing by about thirty per- 
cent from unity for z < 0.8 — 0.9. In this region factoriza- 
tion for these reactions does not hold. On the other hand 
for deep inelastic scattering at Q 2 > 10 GeV 2 the rescat- 
tering of the qq dipole in the photon by the neutron is 
weaker. The different sizes of the absorptive corrections 
in pp compared to j*p interactions prevent a model in- 
dependent extraction of F 2 from the simple factorization 
hypothesis, in agreement with the conclusions reached in 
ref. [tL7|. We also show that the effect of rescattering for 
real photon-proton reactions is important and that one 
can disentangle the difference with highly virtual photons: 
the shift in the neutron energy distributions can be under- 
stood in terms of different absorption effects. 

In addition an extraction of the pion structure func- 
tion is affected by the uncertainty in the determination of 
the pion flux factor. Our fits to the proton induced frag- 
mentation reaction give for the radius parameters Rf c = 
0.2 GeV~ 2 and R c = 0.05 GeV" 2 . More accurate p t neu- 
tron spectra in the case of deep inelastic scattering could 
reduce this uncertainty significantly. At the moment the pt 
spectra in proton induced reactions have been measured 
for too low energies (2a] or too high p t in the relevant z 
range §§. ' 

It may be added that other mechanisms in all three re- 
actions start contributing for z < 0.75; /9-exchangc ||[|,|2l]] 
has been added to the strong proton induced cross section 
presenting a non negligible contribution for 0.5 < z < 0.7. 
In Monte Carlo simulations (LEPTO) the fragmentation 
cross section increases for z < 0.5. Further work is needed 
to understand this part of the cross section theoretically. 

The semi-inclusive reactions in the large z-region show 
a case of transparency reminiscent of diffr active events. 
We think this finding of our paper merits a more accurate 
experimental examination. Exclusive (e.g. p-production ) 
reactions are in seemingly good agreement with theoretical 
calculations based on the dipole picture, which is under- 
lying color transparency. What is important to realize is 
that the final state interactions of virtual partons appear 
correlated with their initial state interaction. 
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Appendix 

Here we give the expression for the function entering the 
Skyrme-type form factor, eq. (H), J38[: 



17. N.N. Niko laev, J. Speth and B.C. Zakharov, hep- 



ph/9708290 



1 + mx m; 



2x — a — b 



b — a 
(39) 

where is the k-th Tchebycheff polynomial of 1st kind, 
with a = 0.003113 b = 280.198488 m = 10~ 8 
a = 6.242336 a x = 4.940353 a 2 = 2.740654 
a 3 = 0.9217577 . 
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